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13.  ABSTRACT  (Maximum  200  Words) 

Neoadjuvant  chemotherapy  is  currently  the  standard  of  care  for  locally  advanced  breast 
cancer  (LABC) .  Monitoring  tumor  response  is  advantageous  for  patients.  This  project  aims 
at  establishing  noninvasive  monitoring  of  neoadjuvant  chemotherapy  in  the  breast  using 
subharmonic  aided  pressure  estimation  (SHAPE)  to  estimate  the  interstitial  fluid  pressure 
(IFP)  in  LABC.  In  vitro  experiments  with  the  ultrasound  contrast  agent  Definity  have 
showed  an  inverse  linear  relationship  between  the  change  in  subharmonic  amplitude  and 
hydrostatic  pressure  (r2=0 . 63-0 . 95,  p<0.01)  over  the  pressure  range  associated  with  breast 
tumors  (0-50  mmHg) .  In  vivo  proof  of  concept  for  SHAPE  as  a  noninvasive  monitor  of  IFP 
(r2=0 . 67-0 . 96,  p<0.01)  has  been  provided  based  on  a  swine  model.  The  work  on  in  vivo  SHAPE 
measurements  of  IFP  was  selected  for  the  final  of  the  American  Institute  of  Ultrasound  in 
Medicine  (AIUM)  2011  Young  Investigator  Award  (for  VG  Halldorsdottir )  and  was  awarded 
FIRST  pi  ace  out  of  90  abstracts  submitted  for  this  competition.  In  total,  97  female,  nude 
rats  were  implanted  with  MDA-MB-231  cells  and  51  animals  were  successfully  studied  with 
SHI  and  SHAPE  21,  24  or  28  days  after  tumor  inoculation.  However,  the  comparisons  between 
SHAPE  and  Stryker  IFP  measurements  have  yet  to  be  completed.  The  SHI  images  of  tumor 
vascularity  were  compared  to  immunohistochemical  markers  of  angiogenesis  (CD31,  COX-2  and 
VEGF)  .  The  strongest  correlation  determined  by  linear  regression  in  this  breast  cancer 
model  was  between  SHI  and  percent  area  stained  with  VEGF  (r=0.53) . 
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4  INTRODUCTION 


In  the  United  States,  close  to  5  -  20  %  of  newly  diagnosed  breast  cancer  and  10  -  30%  of 
all  primary  breast  cancer  is  diagnosed  as  locally  advanced  breast  cancer  (LABC)  [1,  2]. 
Neoadjuvant  chemotherapy  (systemic  preoperative  chemotherapy)  is  currently  the 
standard  of  care  for  LABC  [3,  4].  When  compared  with  adjuvant  chemotherapy 
(postoperative  therapy),  neoadjuvant  chemotherapy  yields  similar  results  for  both  overall 
survival  (70%  for  both)  and  disease-free  survival  (53%  adjuvant,  55%  neoadjuvant)  [5]. 
Thus,  the  postponement  of  surgery  does  not  affect  the  outcome  of  the  treatment  [5,  6],  In 
addition,  neoadjuvant  chemotherapy  offers  considerable  benefits  to  the  patient  as  the 
treatment  can  shrink  the  tumor  and  even  in  some  cases  offer  complete  pathologic 
response  [3,  7].  This  reduction  in  tumor  size  increases  the  possibility  of  breast 
conservation  [3,  5-7].  Maximizing  the  conservation  of  breast  tissue  can  be  of  great 
personal  importance  for  the  self-esteem  and  quality  of  living  of  the  patient  [6]. 
Neoadjuvant  chemotherapy  can  also  offer  an  early  indication  of  the  patient’s  response  to 
chemotherapy.  Consequently,  monitoring  tumor  response  to  neoadjuvant  therapy  gives 
the  possibility  of  adjusting  the  treatment  if  the  patient  is  responding  poorly  or  not  at  all 
resulting  in  substantial  advantages  for  the  patient  [3,  6],  This  project  aims  at  establishing 
noninvasive  monitoring  of  neoadjuvant  chemotherapy  in  the  breast  using  subhannonic 
aided  pressure  estimation  (SHAPE;  U.S.  Patent  6,302,845). 

Generally  interstitial  fluid  pressure  (IFP)  is  10-30  mmHg  higher  in  cancerous  tissue  than 
in  nonnal  tissue  although  values  of  up  to  60  mmHg  have  been  recorded  [8,  9].  Similarly, 
IFP  in  breast  cancer  tumors  has  been  shown  to  be  higher  than  that  of  surrounding  breast 
tissue  [9].  This  increase  is  believed  to  be  due  to  vascularity,  fibrosis  and  difference  in  the 
interstitial  matrix  in  tumors  and  it  can  result  in  poor  transport  of  therapeutic  drugs  to 
tumors  [8].  Taghian  et  al.  used  a  wick-in-needle  technique  to  monitor  the  IFP  of  breast 
cancer  before  and  after  neoadjuvant  chemotherapy  with  two  drugs  used  consecutively 
[10].  When  used  as  a  first  drug  Paclitaxel  decreased  the  IFP  by  36%  (p=0.02)  whereas 
with  Doxorubicin  as  a  first  drug  there  was  only  8%  reduction  (p=0.41).  As  this  was  a 
hypothesis-generating  study  they  did  not  show  any  outcome  related  to  the  relationship 
between  IFP  and  therapy  response  [10].  However,  the  level  of  IFP  has  been  shown  to 
predict  disease  free  survival  for  cervix  cancer  (34%  disease  free  survival  (DFS)  if  IFP  > 
19  mmHg,  68%  DFS  if  IFP  <19  mmHg  (p  =  0.002))  [11],  Thus,  the  level  of  interstitial 
fluid  pressure  (IFP)  in  breast  cancer  tumors  could  potentially  be  used  to  monitor  the 
response  to  neoadjuvant  chemotherapy. 

Contrast  agents  have  been  used  for  two  decades  to  improve  visualization  in  ultrasound 
(US)  imaging  as  they  enhance  the  difference  in  reflectivity  between  tissues  [12],  Because 
of  the  difference  in  compressibility  between  the  medium  and  the  microbubble  any 
changes  in  pressure  induce  changes  in  the  size  of  the  microbubble  [13],  This  in  turn 
affects  the  reflectivity  and  resonant  frequency  of  the  bubble  [13,  14],  In  subharmonic 
imaging  (SHI)  pulses  are  transmitted  at  a  frequency  fo  and  the  echoes  are  received  at  half 
that  frequency  fo/2.  SHI  has  been  showed  to  be  a  feasible  option  for  contrast  enhanced 
imaging  due  to  subhannonic  generation  by  contrast  agents  and  limited  subhannonic 
generation  in  tissues  [15],  Our  group  came  up  with  a  novel  technique,  SHAPE,  utilizing 
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microbubbles  and  the  subharmonic  amplitude  of  the  scattered  signal  [13].  We  showed 
that  there  is  a  linear  relationship  between  the  hydrostatic  pressure  and  the  subhannonic 
amplitude.  We  propose  the  use  of  SHAPE  to  monitor  treatment  response  by 
noninvasively  measuring  the  IFP  in  breast  tumors.  This  offers  several  benefits  to  the 
patient.  As  opposed  to  the  wick-in-needle  technique  SHAPE  is  noninvasive  and  does  not 
inflict  pain.  Furthermore,  it  allows  for  an  early  indication  of  responders  vs.  non¬ 
responders  and  thereby  makes  adjustments  to  therapy  easier.  Moreover,  SHAPE  has  been 
shown  to  have  a  favorable  signal-to-noise  ratio  so  the  subharmonic  amplitude  is  not 
affected  by  background  noise  [13]. 

The  optimal  contrast  agent  and  acoustic  parameters  for  SHAPE  will  be  established  using 
in  vitro  pulse-echo  measurements.  The  SHAPE  algorithm  will  then  be  designed  and 
implemented  on  a  commercial,  state-of-the-art  US  scanner  for  in  vivo  IFP  measurements. 
A  similar  algorithm  has  already  been  set  up  for  cardiac  SHAPE  and  thus  only  a  few 
adjustments  need  to  be  made  to  implement  SHAPE  for  breast  tumors  making  this  very 
cost-effective.  The  in  vivo  experiments  will  be  twofold.  First,  athymic,  nude,  female  rats 
will  be  implanted  with  SKBR3,  MCF-7  or  BT474  human  breast  cancer  cells  and  SHAPE 
used  to  measure  IFP  and  calibrated  by  comparing  the  SHAPE  results  to  IFP 
measurements  obtained  with  an  invasive,  intra-compartmental  pressure  monitor  as  the 
gold  standard.  After  calibration,  human  xenograft  breast  tumors  in  athymic,  nude,  female 
rats  will  be  used  to  evaluate  the  ability  of  SHAPE  to  track  changes  in  IFP  by  studying 
before  and  after  administration  of  a  chemotherapy  agent  (paclitaxel). 

Our  group  has  proposed  that  SHAPE  and  contrast  enhanced  US  imaging  can  be  used  to 
measure  the  IFP  in  LABC  tumors,  thus,  making  it  possible  to  noninvasively  monitor  the 
tumor  response  to  neoadjuvant  chemotherapy.  This  method  would  be  a  considerable 
improvement  from  the  wick-in-needle  technique  currently  used  for  IFP  measurements  in 
LABC  and  allow  for  individualized  treatments  options. 


5  BODY 

The  hypothesis  of  this  project  is  that  IFP  in  breast  tumors  can  be  measured  noninvasively 
using  SHAPE  and  contrast  enhanced  US  thus  improving  the  monitoring  of  neoadjuvant 
chemotherapy.  To  investigate  this  prospect,  in  vitro  pulse-echo  experiments  will  be 
conducted  to  investigate  this  prospect  and  find  the  optimal  contrast  agent  for  SHAPE. 
These  results  will  then  be  used  to  implement  SHAPE  on  a  commercial  scanner.  The 
scanner  will  be  used  for  in  vivo  studies  on  201  rats  with  tumor  xenografts  in  order  to 
calibrate  and  evaluate  SHAPE’S  ability  to  monitor  response  to  neoadjuvant 
chemotherapy.  The  specific  tasks  of  the  project  (as  presented  in  the  original  Statement  of 
Work)  can  be  found  in  Appendix  I. 

First  an  outline  of  the  methods  applied  will  be  given  followed  by  a  presentation  of  the 
results  to  date.  Finally,  the  conclusions  and  future  directions  of  the  research  will  be 
discussed. 
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5. 1  Methods 

In  vitro  experiments 

A  pulse-echo  system,  similar  to  the  one  previously  developed  by  our  group  [13],  was 
constructed  to  test  different  types  of  contrast  agents  for  use  with  high  frequency  SHAPE. 
A  programmable  function  generator  (8116A;  Hewlett  Packard,  Santa  Clara,  CA) 
produced  pulses,  which  were  supplied  to  the  transmit  transducer  after  amplification  in  a 
broadband  50  dB  RF  power  amplifier  (325LA;  ENI,  Rochester  NY).  Received 
microbubble  signals  were  amplified  with  a  low  noise  RF  amplifier  (5052  PR; 
Panametrics,  Waltham,  MA),  digitized  and  processed  (with  the  built-in  FFT  function) 
using  a  digital  oscilloscope  (9350AM;  LeCroy,  Chestnut  Ridge,  NY).  The  amplitude  of 
the  hannonic  and  subharmonic  signal  components  was  obtained  from  spectra  averaged 
over  64  sequences.  Command  delivery  and  data  transfer  was  controlled  by  LabView 
(National  Instruments,  Austin,  TX).  The  subhannonic  amplitude  at  different  static 
pressures  was  measured  using  a  sealed  water  tank  capable  of  withstanding  pressure 
changes  over  200  mmHg.  Single  element  transducers  with  center  frequencies  of  4  to  12 
MHz  were  used  as  transmitter  and  receiver.  The  tank  was  immersed  in  a  water  bath  at 
room  temperature  (25°C).  The  pressure  inside  the  tank  was  monitored  by  a  pressure 
gauge  (OMEGA  Engineering,  Stamford,  CT).  An  inlet  and  an  outlet  on  the  tank  were 
constructed  for  injecting  microbubbles  and  for  applying  extra  hydrostatic  pressure. 

In  order  to  find  the  optimal  contrast  agent  for  high-frequency  SHAPE  changes  in  the 
amplitudes  of  the  first,  second  and  subhannonic  were  measured  for  two  different  contrast 
agents;  Definity  (Lantheus  Medical  Imaging,  N  Billerica,  MA)  and  Sonazoid  (GE 
Healthcare,  Oslo,  Norway)  at  pressures  from  0  to  47  mmHg  (to  simulate  the  IFP  in  breast 
tumors).  Furthermore,  the  frequency  and  acoustic  pressure  were  varied  from  5.7  to  9.3 
MHz  and  0.7- 1.2  MPa,  respectively,  to  detennine  the  optimal  acoustic  parameters.  After 
data  retrieval  the  amplitude  of  the  first,  second  and  subhannonic  was  extracted  using 
MATLAB  (Mathworks,  Natick,  MA).  Three  measurements  were  acquired  at  each  setting 
and  linear  regression  analysis  used  to  detennine  the  relationship  between  hydrostatic 
pressure  and  change  in  amplitude  for  the  first,  second  and  subharmonic.  All  statistical 
analysis  was  conducted  using  Stata  12.0  (Stata  Corporation,  College  Station,  TX). 

Two  alternative  in  vitro  approaches  were  also  investigated.  First,  we  repeated  the  static 
pressure  measurements  using  a  different  -  potentially  more  suitable  -  water  tank.  Small 
(10  ml)  OptiCell  chambers  (Nunc,  Rochester,  NY)  have  been  used  successfully  to 
investigate  interactions  between  contrast  bubbles,  ultrasound  and  cancer  cells  [16]  and 
these  were,  therefore,  selected  as  the  new  container  for  static  pressure  tests.  The  OptiCell 
was  submerged  into  a  larger  water  tank  and  0.2  ml/1  of  Definity  (Lantheus  Medical 
Imaging,  N  Billerica,  MA)  injected.  RF  signals  were  acquired  with  a  Sonix  RP 
ultrasound  scanner  (Ultrasonix,  Richmond,  BC,  Canada)  using  a  high  frequency,  linear 
array  and  compared  to  an  invasive  (needle-in- wick)  pressure  monitor  (Stryker,  Berkshire, 
UK)  as  the  reference  standard.  A  2.5  cm  tissue  mimicking  phantom  was  placed  between 
the  probe  and  the  OptiCell  chamber  to  simulate  tissue  attenuation.  Two  different  transmit 
frequencies  of  6.7  and  10  MHz  were  considered.  Contrast  echoes  were  received  at  half 
the  transmit  frequencies  (i.e.,  3.35  and  5.0  MHz).  The  acoustic  parameters  were  varied 
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and  the  amplitude  of  the  subharmonic  signal  component  was  extracted  and  analyzed  as 
described  above. 

Secondly,  a  completely  new  tank  (with  extra  acoustic  absorbers  incorporated  to  minimize 
reverberations)  was  designed.  The  new  water  tank  was  constructed  and  the  pressure 
inside  the  tank  was  monitored  by  a  pressure  gauge  (OMEGA  Engineering,  Stamford,  CT; 
Fig  1).  A  state-of-the-art  commercial  scanner  Sonix  RP  (Ultrasonix  Medical 
Corporation,  Richmond,  BC,  Canada)  with  a  high  frequency  linear  array  (LI 4-5/3 8)  was 
used  to  acquire  radio-frequency  (RF)  data  at  the  focal  zone  depth  (4.25  cm)  at  a  9  Hz 
framerate  following  injection  of  the  contrast  agent  Definity  (Lantheus  Medical  Imaging, 
N  Billerica,  MA)  in  a  0.2  ml/1  dose.  Two  different  transmit  frequencies  of  6.7  and  10 
MHz  were  considered.  Contrast  echoes  were  received  at  half  the  transmit  frequencies 
(i.e.,  3.35  and  5.0  MHz).  The  acoustic  output  power  was  varied  from  0  to  -20  dB  (in  2 
dB  steps;  0.24-2.05  MPa)  for  a  0  mmHg  hydrostatic  pressure  in  order  to  establish  the 
optimal  sensitivity  for  SHAPE.  Then  the  chamber  pressure  was  varied  from  0  to  50 
mmHg  to  simulate  IFP  in  tumors  and  acoustic  pressure  varied  from  -16  to  -4  dB  (in  4  dB 
steps).  After  RF  data  acquisition,  files  were  transferred  to  a  PC  and  the  amplitude  of  the 
subharmonic  signal  component  was  extracted  using  MATLAB  (Mathworks,  Natick,  MA) 
as  shown  in  Figure  2.  Three  measurements  were  acquired  at  each  setting  and  linear 
regression  analysis  used  to  determine  the  relationship  between  hydrostatic  pressure  and 
change  in  subharmonic  amplitude.  All  statistical  analyses  were  conducted  using  Stata 
9.0-12.0  (Stata  Corporation,  College  Station,  TX). 

Moreover,  a  novel,  simulation  model  of  the  dynamics  of  an  encapsulated  microbubble 
contrast  agent,  developed  as  part  of  a  previous  DOD  supported  project  [14],  was 
modified  in  order  to  account  for  ambient  pressure  variations  and  different  shell 
parameters  to  establish  the  optimal  contrast  microbubble  for  SHAPE.  A  nonlinear 
extension  of  the  original  viscoelastic  model  was  pursued  by  considering  a  quadratic 
elasticity  model  where  the  interfacial  elasticity  vary  linearly  with  area  fraction  as  well  as 
an  exponential  model  (i.e.,  the  elasticity  vary  exponentially). 


Figure  1.  The  re-designed  water  tank  scanned  by  the  L14-5  transducer.  Notice  the 
digital  pressure  gauge  on  the  top  of  the  tank. 
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Figure  2.  Block  diagram  of  the  off-line  processing  performed  on  the  RF  data. 

In  vivo  experiments 

All  research  studies  involving  animals  were  approved  by  the  Institutional  Animal  Care 
and  Use  Committee  of  Thomas  Jefferson  University  and  conducted  in  accordance  with 
the  guidelines  provided  by  the  National  Institutes  of  Health. 

Our  group  has  worked  in  partnership  with  Ultrasonix  Medical  Corporation  to  implement 
SHAPE  for  cardiac  use  on  a  state-of-the-art  commercial  scanner  Sonix  RP  (Ultrasonix 
Medical  Corporation,  Richmond,  BC,  Canada)  with  a  phased  array  (PA4-2).  Several 
experiments  have  been  carried  out  in  canines  to  investigate  cardiac  SHAPE  supported  by 
funding  from  the  AHA  and  the  NIH.  RF  data  from  these  experiments  was  analyzed  off¬ 
line  using  Matlab.  The  software  developed  by  our  group  for  this  analysis  is  not  site- 
specific  and  will  also  be  used  to  analyze  the  data  we  will  acquire  from  in  vivo  breast 
SHAPE  in  rats  as  part  of  this  project. 

Furthermore,  an  opportunity  to  provide  a  proof-of-concept  of  the  use  of  SHAPE  for 
estimating  IFP  and  test  the  invasive  (needle  based)  Stryker  pressure  monitoring  system 
(the  reference  standard)  presented  itself.  As  part  of  an  ongoing  NIH  study,  a  unique, 
naturally  occurring  tumor  model,  the  Sinclair  swine  with  melanoma,  was  being  studied. 
We  obtained  IFP  measurements  from  the  tumors  and  surrounding  tissue  using  the  Stryker 
needle  based  system,  which  provided  us  with  a  chance  to  assess  the  dependence  of  this 
technique  on  the  angle  between  the  needle  and  the  tissue.  Moreover,  subharmonic 
signals  were  acquired  during  an  infusion  of  Definity  (6.25  ml/l/min)  with  the  Sonix  RP 
and  a  linear  array.  Data  were  obtained  at  6.7  and  10  MHz  (i.e.,  subharmonic  frequencies 
of  3.75  and  5.0  MHz,  respectively)  with  acoustic  outputs  of  -4,  -8  and  -12  dB  and 
analyzed  further  during  the  last  year.  Five  (5)  swine  were  studied  (one  melanoma  per 
swine)  at  no  cost  to  this  project. 

Originally,  we  planned  to  study  three  human  breast  cancer  cell  lines  SKBR3,  MCF-7  and 
BT474  in  this  project.  However,  our  collaborators  advised  us  against  trying  to  work  with 
so  many  different  cell  lines  in  the  short  time  available  for  the  in  vivo  rat  experiments. 
Hence,  it  was  decided  to  focus  our  efforts  on  the  MDA-MB-231  breast  cancer  cell  line, 
since  it  is  known  to  grow  successfully  in  rats  [17-19].  This  cell  line  was  also  selected, 
because  it  is  able  to  upregulate  intracellular  pH  to  near  normal  levels  and,  therefore,  have 
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a  higher  likelihood  of  remaining  viable  under  the  adverse  (low  pH)  conditions  that  are 
typical  of  the  tumor  microenvironment  [20].  Cells  in  such  an  environment  are  the  most 
stimulatory  in  terms  of  inducing  tumor  angiogenesis.  The  MDA-MB-23 1  cell  line  was 
cultured  in  Eagle’s  Minimum  Essential  medium  alpha  modification  (alpha  MEM;  Sigma- 
Aldrich  Co.,  St.  Louis,  MO)  in  an  incubator  at  37°  C  with  a  5  %  COo  in  air  atmosphere. 

Female,  5-8  weeks  old,  nude,  athymic  rats  (Charles  River  Laboratories,  Wilmington, 
MA),  weighing  between  150  and  250  g,  were  injected  subcutaneously  with 
approximately  5  x  106  tumor  cells  into  the  right  mammary  fat  pad.  However,  it  took 
almost  8  months  to  arrive  at  a  stable  and  consistent  level  of  cell  growth  for  the  MDA- 
MB-23 1  cells,  which  also  produced  tumor  growth  in  vivo.  Twelve  (12)  rats  were 
evaluated  in  this  period  and  our  studies  indicated  that  time  points  of  21,  24  or  28  days 
after  tumor  inoculation  ensured  enough  difference  in  tumor  volume  and  tumor 
vasculature  for  effectively  monitoring  the  progression  of  angiogenesis 

For  the  SHAPE  experiments,  it  is  critical  to  select  the  appropriate  acoustic  output  (AO) 
power  (i.e.,  in  the  growth  phase  of  the  subharmonic  signal  generation).  Hence,  our  group 
developed  an  automated  AO  optimization  algorithm  to  identify  the  AO  level  eliciting 
growth  stage  subharmonic  signals  to  help  eliminate  the  problems  of  acquiring  and 
analyzing  the  data  offline  at  all  AO  levels  (as  in  our  previous  studies)  [21].  The  AO 
optimization  algorithm  was  implemented  on  a  Sonix  RP  scanner  along  with  pulse 
inversion  SHI  transmitting  at  8  MHz  and  receiving  at  4  MHz.  The  optimization 
algorithm  stepped  the  ultrasound  scanner  from  0  to  100  %  AO.  A  logistic  equation  fitting 
function  was  applied  with  the  criterion  of  minimum  least  squared  error  between  the  fitted 
subharmonic  amplitudes  and  the  measured  subharmonic  amplitudes  as  a  function  of  the 
AO  levels  and  the  optimum  AO  level  was  chosen  corresponding  to  the  inflection  point 
calculated  from  the  fitted  data  (Fig.  3). 


Acoustic  Power  vs.  Subharmonic  Signal 


Figure  3.  Example  of  the  outcome  of  the  AO  optimization  algorithm. 


9 


For  the  initial  calibration  studies,  25  rats  received  MDA-MB-23 1  cell  line  implantations. 
The  IFP  of  the  nonnal  mammary  fat  pad  as  well  as  the  tumor  was  measured  with  an  intra- 
compartmental  pressure  monitor  (selected  as  the  reference  standard  technique  for  this 
study).  The  Stryker  pressure  monitor  is  an  accurate  but  invasive  technique  for  sampling 
the  interstitial  fluid  and  measuring  the  IFP.  Grayscale  SHI  was  performed  with  the 
modified  Sonix  RP  scanner  using  the  LI 4-5  linear  array  after  contrast  administration 
(Definity;  dose:  36  pi)  and  the  largest  diameter  of  the  tumor  was  imaged.  Maximum 
intensity  projections  of  3  SHI  measurements  were  constructed  off-line  and  averaged  for 
each  SHAPE  data  point.  Furthermore,  as  direct  measurement  of  the  tumor  IFP  using  a 
needle  based  system  will  vary  due  to  cancer  heterogeneity,  all  tumor  IFP  data  were  also 
averaged  over  3  measurements. 

The  remaining  60  nude  rat  xenografts  were  used  to  establish  the  ability  of  SHAPE  to 
monitor  therapy  responses.  SHAPE  and  the  Stryker  pressure  monitor  measurements  were 
perfonned  before  and  after  administration  of  a  single  IV  injection  of  5  mg/kg  of  the 
chemotherapy  agent  paclitaxel  (Mayne  Pharma,  Paramus,  NJ)  at  the  selected  post¬ 
implantation  time  points.  This  is  not  similar  to  the  therapy  regimes  employed  in  humans 
(which  typically  involve  a  40  mg/kg  dose),  but  the  literature  clearly  describes  how  the 
human  dose  (even  weight-based)  is  unsuitable  for  rats  [22].  Hence,  the  choice  of  5 
mg/kg  as  the  treatment  dose  utilized.  Paclitaxel  is  known  to  reduce  tumor  IFP.  SHAPE 
measurements  were  perfonned  pre  and  48  hours  post  therapy,  since  this  is  when  the 
biggest  reduction  in  IFP  is  expected.  Additionally,  tumor  volumes  were  measured  on 
days  21,  24  and  28  to  allow  growth  curves  to  be  established  and  compared.  The  entire 
study  lasted  no  longer  than  1  hour. 

At  the  end  of  the  US  study,  rats  were  euthanized  by  placing  them  in  a  gas  chamber  and 
saturating  the  air  with  CCL.  Tumors  were  then  surgically  removed  and  scanned  ex  vivo 
to  identify  the  US  imaging  plane  studied  in  vivo.  Careful  attention  was  paid  to  the 
labeling  of  each  specimen  by  specimen-ink  marking  the  true  front,  back,  right  and  left 
[23].  Then  specimens  were  placed  in  10  %  formalin  phosphate  (Fisher  Scientific, 
Houston,  TX)  for  12  to  24  hours  to  fix  all  angiogenic  markers,  dehydrated  in  graded 
alcohols,  cleaned  in  xylene,  and  embedded  in  paraffin  using  standard  methods.  The 
process  of  angiogenesis  involves  the  activation,  migration,  and  proliferation  of 
endothelial  cells  [24-25].  Vascular  endothelial  growth  factor  (VEGF)  is  regarded  as  a 
major  regulator  of  tumor-associated  angiogenesis,  since  it  promotes  tumor  growth, 
invasion,  and  metastasis  [24-27].  Cyclooxygenase-2  (COX-2)  is  another  positive 
regulator  of  angiogenesis,  which  is  expressed  in  response  to  growth  factors,  tumor 
promoters,  or  cytokines  [28-29].  Platelet  endothelial  cell  adhesion  molecule  (PEC AM  or 
CD31)  is  a  potent  endothelial  cell  marker,  since  its  expression  is  restricted  to  vascular 
system  platelets  and  endothelium  [30],  Each  specimen  was  immunohistochemical 
stained  using  a  monoclonal  antibody  against  PECAM  (anti-CD31;  Dako  Corporation, 
Carpinteria,  CA),  a  polyclonal  antibody  against  COX-2  (Santa  Cruz  Biotechnology, 
Santa  Cruz,  CA),  and  a  monoclonal  antibody  against  VEGF  (Oncogene  Research 
Products,  San  Diego,  CA)  as  reference  standards.  Finally,  the  sections  were  mounted 
onto  glass  slides  for  analysis. 
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A  semi-automated  histomorphometry  system  incorporating  a  DXC-970MD  color  CCD 
camera  (Sony  Corporation,  Tokyo,  Japan)  connected  to  an  SMZ-10A  microscope  (Nikon, 
Melville,  NJ,  USA)  was  used  to  analyze  both  SHI  images  and  specimens.  The 
histomorphometry  system  employed  a  4x  objective  and  lOx  ocular  magnification  (total 
magnification:  40x)  to  provide  a  digital  view  of  each  histological  slide  on  a  desktop 
computer  [23],  Area  measurements  of  stains  and  enhanced  pixels  relative  to  the  total 
tumor  area  were  obtained  using  ImagePro  Plus  software  (Media  Cybernetics,  Silver 
Spring,  MD).  Briefly,  the  tumor  was  manually  outlined  on  the  images  and  the  RGB 
channels  used  to  automatically  segment  and  count  the  number  of  enhanced  or  stained 
pixels  as  well  as  the  total  number  of  pixels  within  the  tumor.  Fractional  tumor 
vascularity  (FV)  was  calculated  (in  percent)  as  the  number  of  contrast  enhanced  pixels 
(Cj)  relative  to  the  total  number  of  pixels  (c,  +  x,): 

Yc,. 

FV  =  x  100%  m 

Results  were  compared  using  linear  regression  analysis. 


5.2  Results  and  Discussion 

In  vitro  experiments 

Over  the  pressure  range  of  0  -  50  mmHg  (simulating  the  IFP  in  breast  tumors)  both 
Definity  and  Sonazoid  showed  an  inverse  linear  relationship  between  the  change  in 
subhannonic  amplitude  and  hydrostatic  pressure  (r  =  0.76-0.91  ,  P  <  0  .01),  when  studied 
with  the  pulse-echo  setup.  These  results  are  consistent  with  previous  results  reported  by 
our  group  [13].  An  example  of  the  subhannonic  amplitude  for  Sonazoid  at  0  mmHg  and 
47  mmHg  can  be  seen  in  Figure  4  (transmitting  frequency  7.5  MHz  and  acoustic  pressure 
0.7  MPa).  This  process  was  markedly  more  time-consuming  then  originally  envisaged, 
due  to  difficulty  in  aligning  the  single  element  transducers  at  these  higher  frequencies 
(compared  to  our  original  efforts  e.g.,  in  [13])  and  this  issue  delayed  the  project  by 
approximately  6  months. 

As  an  alternative  OptiCell  measurements  with  Definity  were  conducted  and  these  showed 
an  inverse  linear  relationship  between  the  change  in  subhannonic  amplitude  and 
hydrostatic  pressure  (r2  =  0.79-0.99,  p  <  0.01).  However,  the  decrease  in  subhannonic 
amplitudes  recorded  for  pressure  changes  from  0  to  50  mmHg  varied  from  13.7  to  22.9 
dB  and  16.4  to  21.8  dB  (depending  on  acoustic  pressures)  for  6.7  and  10  MHz 
transmission  frequencies,  respectively.  These  levels  of  subhannonic  amplitude  variation 
were  markedly  higher  than  our  previous  results  obtained  over  a  much  larger  pressure 
range  (up  to  200  mmHg)  [13,  31-32].  After  some  extensive  testing,  we  discovered  that 
standing  waves  markedly  influenced  results  and  made  measurements  impossible  to 
reproduce.  Hence,  our  use  of  the  OptiCell  setup  had  to  be  abandoned  and  a  completely 
new  tank  (with  extra  acoustic  absorbers  incorporated)  had  to  be  designed  and  built  (as 
described  in  Section  5.1).  This  issue  delayed  the  project  further  and  we,  therefore, 
requested  and  were  granted  a  one  year  no  cost  extension. 
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Figure  4.  Comparison  of  the  subharmonic  amplitude  (circle)  measured  with  Sonazoid  at 
a)  0  mmHg  and  b)  47  mmHg pressures. 


Once  the  new  water  tank  (Fig.  1)  was  available,  we  achieved  much  more  reproducible 
results.  The  characteristic  sigmoidal  curve  for  the  relationship  between  the  subharmonic 
amplitude  and  the  acoustic  power  showed  the  three  stages  of  subhannonic  generation: 
i.e.,  occurrence  (-20  dB  to  -16  dB)  where  there  was  minimal  change  in  the  subharmonic 
amplitude;  growth  (-16  dB  to  -  4  dB)  where  there  was  a  sharp  rise  in  subharmonic 
amplitude  and  the  subharmonic  is  most  sensitive  to  pressure  changes,  and  finally 
saturation  (-4  dB  to  0  dB)  where  again  there  is  minimal  subharmonic  amplitude  change 
(Fig.  5).  Based  on  this  relationship  four  acoustic  power  levels  covering  the  growth  stage 
and  its  boundaries  were  chosen  for  further  investigation  of  hydrostatic  pressure 
estimation  using  SHAPE;  -16,  -12,  -8  and  -4  dB.  These  values  correspond  to  0.33,  1.06, 
1.33,  1.68  MPa  and  0.24,  1.21,  1.52  and  1.78  MPa  peak  to  peak  for  6.7  MHz  and  10 
MHz,  respectively,  as  measured  with  a  0.2  mm  needle  hydrophone  (Precision  Acoustics, 
Dorchester,  Dorset,  UK). 


In  the  water-tank  an  inverse  relationship  was  seen  over  a  pressure  change  from  0  to  50 
mmHg  at  both  6.7  MHz  and  10  MHz  transmission  frequencies  and  all  acoustic  power 
levels  -16  dB  to  -4  dB  (r2  >  0.63;  p  <  0.05;  Fig.  6).  The  r2  values  from  the  linear 
regression  analysis  were  consistently  higher  for  a  transmission  of  10  MHz  (receiving  the 
subharmonic  at  5  MHz)  than  for  6.7  MHz  (receiving  at  3.35  MHz).  Moreover,  pressure 
estimates  obtained  at  -16  dB  (start  of  the  growth  phase)  showed  limited  sensitivity 
compared  to  the  higher  acoustic  power  settings  likely  due  to  a  lack  of  subharmonic 
generation  (still  in  the  occurrence  phase).  The  largest  drop  in  subharmonic  amplitude 
(corresponding  to  the  maximum  sensitivity  for  pressure  estimation),  11.36  dB  over  50 
mmHg,  was  seen  at  10  MHz  and  -8  dB  (r“  =  0.95;  p  <  0.01;  Fig.  7).  This  is  consistent 
with  previous  results  reported  by  our  group  [13]  and  these  efforts  represent  the 
fulfillment  of  tasks  la,  Id,  If,  2a  and  2c  in  the  original  Statement  of  Work  (SOW). 
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6.7  MHz 
-10  MHz 


Figure  5.  Subharmonic  response  to  changes  in  acoustic  power  with  the  occurrence, 
growth  and  saturation  phases  clearly  seen. 
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r2:  0.63-0.90 
p  <0.05 

■  10  MHz 

r2:  0.89-0.95 
p  <0.05 


Figure  6.  Maximum  decrease  in  subharmonic  signal  amplitude  for  Definity  as  a 
function  of  frequency  and  acoustic  power  (N  =  3). 


The  previously  developed  simulation  model  [14]  was  modified  to  include  nonlinear 
extensions  of  the  viscoelasticity  and  this  model  has  now  been  published  [33].  The  intent 
was  to  better  account  for  the  experimentally  observed  changes  in  subharmonic  signal 
amplitudes  as  a  function  of  hydrostatic  pressures.  The  model  showed  that  the 
determining  parameter  of  subharmonic  response  is  the  ratio  of  the  excitation  to  the 
resonance  frequency.  Changing  the  ambient  pressure  changes  the  resonance  frequency 
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Figure  7.  Inverse  linear  relationship  between  pressure  and  subharmonic 
signal  amplitude  for  Definitv  at  10  MHz  and  -8  dB  acoustic  output  power 
(N  =  3). 


and  thereby  the  frequency  ratio.  For  different  acoustic  excitation  pressure  levels, 
changing  ambient  pressure  can  either  increase  or  decrease  the  subharmonic  response 
depending  on  this  ratio.  For  some  range  of  parameters,  the  variation  is  far  more 
complicated.  This  behavior  is  clearly  at  odds  with  the  experimental  observations 
mentioned  above  and  is  currently  not  well  understood.  These  discrepancies  may  be  due 
to  encapsulation  buckling  or  strain  softening  in  describing  nonlinear  oscillations, 
specifically  the  subharmonic  response  [33].  This  investigation  has  been  published  [33- 
34].  This  constitutes  the  completion  of  tasks  lb,  lc  and  le. 

In  vivo  experiments 

Software  to  analyze  RF  data  from  the  Sonix  RP  scanner  has  been  improved  and  the  best 
method  to  extract  the  subharmonic  signal  components  from  the  frequency  spectrum  has 
been  established.  This  software  is  applicable  to  RF  data  acquired  from  any  tissues 
(including  breast  tumors).  However,  due  to  the  availability  of  canine  cardiac  data 
(obtained  with  funds  from  other  sources,  but  processed  with  funds  from  this  grant),  our 
initial  work  focused  on  cardiac  SHAPE  pressure  estimates  [32,  35-37].  Briefly,  the 
unprocessed  RF  data  for  each  accumulated  pulse  (Fig.  8A)  was  transformed  to  the 
Fourier  domain  and  the  subharmonic  signal  amplitude  (at  half  the  fundamental  frequency 
-  Fig.  8B)  was  extracted  as  the  average  signal  in  a  40  %  bandwidth  around  the 
subharmonic  frequency  (i.e.,  here  1.25  MHz).  The  extracted  subharmonic  signals  from 
all  pulses  were  processed  using  a  moving  average  filter  to  eliminate  noise  spikes.  The 
range  of  the  subharmonic  signal  (i.e.,  the  difference  between  maximum  and  minimum 
subharmonic  amplitude)  was  compared  from  each  pulse  contour  (after  eliminating 
excessively  noisy  pulses)  for  each  incident  acoustic  pressure.  Since  the  pressure  tracking 
using  microbubbles  is  an  incident  acoustic  pressure  dependent  phenomenon,  the  incident 
acoustic  pressure  with  maximum  stable  subharmonic  range  was  then  selected  for  LV 
pressure  tracking;  as  shown  in  Figure  8C. 
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Figure  8:  Steps  involved  in  extracting  and  processing  the  subharmonic  signal  for 
SHAPE.  A  typical  signal  obtained  with  pulsed  Doppler  (A)  and  the  frequency  domain 
representation  (B)  of  the  signal  in  (A).  The  processed  subharmonic  signal  from  all  the 
pulses  (solid  line)  and  the  pressure  catheter  data  (dotted  line)  are  shown  in  (C).  Note  the 
inverse  relationship  between  the  subharmonic  signal  and  the  pressure  obtained  via  the 
pressure  catheter  (in  agreement  with  our  previously  published  in  vitro  results). 

Results  from  seven  canines  indicate  that  maximum  absolute  errors  on  the  order  of  -1.36 
to  2.50  mmHg  can  be  obtained  for  diastolic  left  ventricular  pressures  if  the  aortic  pulse 
pressure  values  are  known;  otherwise  maximum  absolute  errors  for  LV  estimates 
compared  to  the  reference  standard  were  slightly  larger  (up  to  9  mmHg)  [35-37]. 
Moreover,  estimating  the  peak  systolic  RV  pressure  resulted  in  errors  of  0.0  to  3.5  mmHg 
[35-37].  Moreover,  a  patent  application  has  been  submitted  to  the  U.S  Patent  Office 
based  on  these  development  efforts  [38],  This  effort  represents  the  conclusion  of  task  2b. 

A  unique  opportunity  to  test  the  Stryker  pressure  monitoring  system  and  provide  proof- 
of-concept  of  the  use  of  SHAPE  for  estimating  IFP  was  pursued  in  the  Sinclair  swine 
melanoma  model.  Five  (5)  swine  were  studied,  but  one  of  the  swine  was  eliminated  from 
the  study,  due  to  technical  difficulties.  The  difference  between  the  subharmonics  and  the 
IFP  in  tumors  and  their  surrounding  tissues  was  most  linear  for  10  MHz  transmission  (r2 
>  0.67;  p  <  0.05;  Fig.  9)  [39-40].  The  difference  between  tissue  and  tumor  IFP  is  rather 
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large  and  may  be  caused  by  the  sensitivity  of  the  Stryker  system  to  the  angle  between  the 
needle  and  the  tissue  being  studied  (this  issue  is  currently  being  investigated  further). 
The  slope  from  3  out  of  4  animals  was  very  similar  to  each  other  (-0.19  ±  0.07)  and  to  the 
in  vitro  slope  (-0.22),  while  the  last  swine  (i.e.,  swine  1)  showed  a  large  spread  within  the 
normal  tissue  IFP.  In  conclusion,  these  results  provide  proof  of  concept  for  the  feasibility 
of  using  SHAPE  as  a  noninvasive  monitor  of  IFP. 

Finally,  the  work  on  in  vivo  SHAPE  measurements  of  IFP  in  melanoma  swine  [39]  was 
selected  for  the  final  of  the  American  Institute  of  Ultrasound  in  Medicine  (AIUM)  2011 
Young  Investigator  Award  (for  V.  G.  Halldorsdottir)  and  was  awarded  FIRST  place  out 
of  90  abstracts  submitted  for  this  competition. 

Of  the  25  rats  implanted  with  MDA-MB-231  for  the  initial  calaibration  studies,  16  (64 
%)  exhibited  tumor  growth  and  13  were  successfully  imaged  in  completion  of  task  3  a  and 
3b.  In  the  treatment  groups  (for  days  21,  24  and  28)  60  rats  received  injections  of  the 
breast  cancer  cells  into  the  right  mammary  fatpad,  which  resulted  in  38  animals  with 
tumors  and  amongst  those  26  were  successfully  imaged.  The  15  rats  with  tumor  growth 
that  failed  the  imaging  study  were  mainly  for  technical  reasons  (the  experimental 
software  on  the  Sonix  RP  crashed  occasionally).  The  treatment  groups  therefore  ended 
up  containing  8,  9  and  9  rats  for  21,  24  and  28  days  post-implantation,  respectively.  In 
each  group  there  were  2  rats  that  were  designated  as  control  cases  (i.e.,  they  did  not 
receive  injections  of  paclitaxel).  SHI  depicted  the  tortuous  morphology  of  tumor 
neovessels  and  delineated  areas  of  necrosis  better  than  fundamental  ultrasound  imaging, 
due  to  the  suppression  of  tissue  signals.  However,  the  degree  of  tissue  suppression  was 
less  than  other  in  vivo  studies  undertaken  by  our  group  [31,  41-42],  which  was  attributed 
to  the  nonlinear  transfer  function  of  the  scanner. 


Interstitial  Fluid  Pressure  [mmHg] 


♦  Swine  1:  -4  dB 
■  Swine  2:  -8  dB 
A  Swine  3:  -8  dB 

•  Swine  4:  -4  dB 


70 


Figure  9.  In  vivo  measurements  showing  SHAPE  results  compared  to 
the  pressure  monitor.  The  relative  difference  between  tissue  and  tumor 
IFP  is  clearly  captured  by  SHAPE. 
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The  change  in  tumor  volume  as  well  as  in  IFP  (measured  with  the  Stryker  reference 
standard)  from  pre  to  post  treatment  is  depicted  in  Figure  10.  The  change  encountered 
for  the  control  animals  have  also  been  included.  For  both  markers  the  standard 
deviations  are  quite  large,  but  for  the  IFP  values  a  clear  trend  is  seen  with  the  IFP 
reducing  more  as  time  increases.  However,  for  both  tumor  volume  and  IFP 
measurements  there  were  no  statistically  significant  differences  between  the  control  and 
the  treatment  groups  (p  >  0.37).  It  should  be  pointed  out  that  the  main  analysis  planned 
for  this  project  -  the  comparisons  between  SHAPE  and  Stryker  IFP  measurements  -  is 
still  ongoing.  This  is  due  to  both  the  students  on  this  project  (V  Halldorsdottir  and  A 
Marshall)  becoming  parents  within  the  last  4  months  and  coupled  with  the  issues  of 
pregnancy,  they  were  not  able  to  finalize  the  data  analysis  in  time.  This  work  (under  task 
3d)  is  still  ongoing. 

The  SHI  images  were  rendered  into  maximum  intensity  projection  (MIP)  mode  in  order 
to  better  isolate  contrast  enhanced  pixels  from  tissue  background  and  noise  signals  (Fig. 
11)  when  calculating  the  FV.  When  analyzing  the  immunohistochemical  markers  of 
angiogenesis,  VEGF  was  found  to  vary  significantly  over  time  (p<0.031),  while  the 
strongest  correlation  determined  by  linear  regression  in  this  breast  cancer  model  was 
between  SHI  and  percent  area  stained  with  VEGF  (r  =  0.53;  Table  1).  This  represents  the 
completion  of  Task  3c.  It  should  be  pointed  out  that  the  results  in  Table  1  are  based  on 
32  rats  and  involve  values  obtained  over  the  entire  tumor  area.  However,  given  the 
known  tumor  heterogeneity,  it  is  very  likely  the  averaging  over  the  entire  area  will  not  be 
optimal  [23].  Hence,  we  intend  to  re-analyze  this  data  split  into  smaller  regions  of 
interest.  Another  limitation  is  the  inherent  problem  in  matching  ultrasound  images  with 
beam-widths  on  the  order  of  mm  to  pathology  specimens  cut  in  5  pm  sections;  even 
though  great  care  was  taken  to  match  the  imaging  planes  to  the  pathology  sections  (as 
described  above). 

Finally,  the  work  on  comparing  in  vivo  SHI  measurements  of  tumor  vascularity  in  murine 
xenografts  to  immunohistochemical  markers  [43]  was  selected  for  the  final  of  the  AIUM 
2013  Young  Investigator  Award  (for  A.  Marshall).  Although  placements  beyond  first 
place  are  not  released,  it  should  be  noted  that  this  abstract  was  one  out  of  1 1  selected 
from  a  pool  of  1 1 1  abstracts  submitted  for  this  competition. 
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Figure  10.  The  change  in  tumor  volume  (A)  and  IFP  measurements  (B)  from  pre  to  post 
administration  of  Paclitaxel.  Notice  how  the  treatment  results  in  lower  IFP  values  at 
days  24  and  28,  as  expected. 
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Figure  11.  Breast  tumor  xenograft  (arrows)  depicted  in  maximum  projection  intensity 
(MIP)  SHI  mode  (A)  and  the  pixels  corresponding  to  vascularity  in  red  within  the  tumor 
boundary  (green  dotted  line)  shown  in  (B). 


Table  1.  Correlation  between  SHI  and  angiogenic  markers  in  32  rats. 


VEGF  vs.  SHI 

r  p 

COX-2  vs.  SHI 

r  p 

CD-31  vs.  SHI 

r  p 

Day  21  (n=ll) 

0.202 

0.551 

0.158 

0.642 

0.304 

0.363 

Day  24  (n=ll) 

0.529 

0.094 

0.134 

0.692 

0.114 

0.737 

Day  28  (n=10) 

0.032 

0.943 

0.063 

0.857 

0.045 

0.901 

Overall  (n=32) 

0.032 

0.890 

0.071 

0.689 

0.071 

0.689 

6  KEY  RESEARCH  ACCOMPLISHMENTS 

•  SHAPE  experiments  were  conducted  with  Definity  and  Sonazoid  in  vitro  in  a 
pressurized  watertank  (range:  0-50  mmHg)  using  the  pulse-echo  setup. 

•  SHAPE  experiments  were  conducted  in  vitro  with  Definity  using  the  Sonix  RP 
scanner  and  two  different  probes;  LI 4-5/3 8  and  L9-4/38. 

•  Hydrostatic  pressure  is  inversely  related  to  the  change  in  subharmonic  amplitude 
(r2  =  0.76-0.91,  p<  0.01). 
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•  A  computer  model  to  simulate  the  behavior  of  microbubbles  as  a  function  of 
pressure  has  been  developed. 

•  SHAPE  experiments  were  conducted  in  vitro  in  a  new,  improved  water-tank. 

•  The  maximum  decrease  obtained  with  Definity  over  50  mmHg  was  11.36  dB 
(r2=0.95;  p<0.01)  at  10  MHz  and  -8  dB  power. 

•  Software  for  processing  of  in  vivo  SHAPE  data  has  been  developed  and 
optimized. 

•  In  vivo  proof  of  concept  for  SHAPE  as  a  noninvasive  monitor  of  IFP  has  been 
provided  (r2  =  0.67  -  0.96,  p  <  0.01). 

•  The  L9-4  linear  array  has  been  selected  for  the  in  vivo  studies. 

•  Sonix  RP  scanner  design  requirements  have  been  developed  and  implemented. 

•  Lantheus  Medical  decided  to  support  this  project  with  230  vials  of  Definity  for 
use  in  the  in  vivo  experiments. 

•  In  vivo  SHAPE  data  from  26  rats  before  and  after  treatment  with  paclitaxel  have 
been  acquired,  but  not  yet  analyzed. 

•  Initial  data  analysis  indicate  that  SHI  measures  of  vascularity  may  correlate  with 
the  expression  of  VEGF  (r  =  0.53). 
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Thirty-Sixth  International  Symp.  on  Ultrasonic  Imaging  and 
Tissue  Characterization,  Washington  DC,  USA. 

•  On  the  utility  of  subhannonic  contrast  microbubble  signals. 

•  Dynamics  of  contrast  microbubbles  and  their  subhannonic 

response  for  noninvasive  blood  pressure  estimation 

Era  of  Hope,  Dept  of  Defense  Breast  Cancer  Research 
Meeting,  Orlando,  FL,  USA. 

•  Subhannonic  aided  pressure  estimation  for  monitoring 
neoadjuvant  chemotherapy  of  locally  advanced  breast 
cancer  (poster). 

IEEE  2011  International  Ultrasonics  Symposium,  Orlando,  FL, 
USA. 
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•  Dual  grayscale  and  subharmonic  ultrasound  imaging  on 
modified  ultrasound  scanner  in  2  and  3D  (poster). 

November  27  -  The  97th  Scientific  Assembly  and  Annual  Meeting  of  the 

December  2,  201 1  Radiological  Society  of  North  America,  Chicago,  IL,  USA. 

•  Subhannonic  aided  pressure  estimation  for  monitoring 
interstitial  fluid  pressure  in  swine  melanomas:  initial  in 
vitro  and  in  vivo  results. 

•  Dual  grayscale  and  subharmonic  US  imaging  on  a  modified 
commercial  scanner  (eposter). 

•  Processing  acoustic  data  from  US  contrast  agents  for 
ambient  pressure  estimation  (eposter). 

•  Noninvasive  cardiac  pressure  estimation  using  subhannonic 
microbubble  signals  (eposter). 

November  29,  2011  QIBA  Ultrasound  Quantitative  Biomarker  Exploratory 

Meeting,  Chicago,  IL,  USA. 

•  Pressure  measurements  with  contrast  agents. 

February  3,  2012  Reunion  Intemacional  Sobre  Avances  en  Hepatologia  IV, 

Valladolid,  Spain. 

•  Subhannonic  aided  pressure  estimation  (SHAPE):  a 

noninvasive  technique  for  the  measurement  of  fluid 
pressures  using  ultrasound  contrast  agents. 

October  7  -  10,  2012  IEEE  2012  International  Ultrasonics  Symposium,  Dresden, 

Germany. 

•  Noninvasive  right  ventricular  pressure  estimation  in  vivo 
using  the  subhannonic  emissions  from  ultrasound  contrast 
agents. 

April  6  -  10,  2013  The  58th  Annual  Convention  of  the  American  Institute  of 

Ultrasound  in  Medicine,  New  York  City,  NY,  USA. 

•  Subhannonic  imaging  of  angiogenesis  in  a  murine  breast 

cancer  model 


V.  G.  Halldorsdottir  was  selected  as  a  finalist  for  the  AIUM  2011  Young  Investigator 
Award  based  on  her  work  in  the  abstract  labled  *  above.  The  competition  was  judged  at 
the  Annual  Conference  of  the  AIUM  (in  April  2011  in  New  York  City,  NY)  and  was 


awarded  FIRST  place  out  of  90  abstracts  submitted  for  this  competition. 


A.  Marshall  was  selected  as  a  finalist  for  the  AIUM  2013  Young  Investigator  Award 
based  on  his  work  in  the  abstract  labled  **  above.  The  competition  was  judged  at  the 
Annual  Conference  of  the  AIUM  (in  April  2013  in  New  York  City,  NY)  and  although  it 
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Patents 

F.  Forsberg,  J.  K.  Dave,  V.  G.  Halldorsdottir,  J.  R.  Eisenbrey.  Methods  for  improved 
selection,  processing  and  display  of  subharmonic  microbubble  signals  as  pressure 
estimates.  U.S.  Patent  pending  61/498,278,  June,  2011. 

Degrees 

While  no  students  have  yet  graduated  based  on  the  data  acquited  as  part  of  this  project,  our 
expectation  currently  is  that  one  (1)  PhD  and  two  (2)  MSc’s  will  be  granted  by  Drexel 
University. 


8  CONCLUSIONS 

We  have  investigated  six  US  contrast  agents  for  use  in  SHAPE  and  shown  that  Dcfinity 
(Lantheus  Medical  Imaging,  N  Billerica,  MA)  has  an  inverse  linear  relationship  between 
the  change  in  subharmonic  amplitude  and  hydrostatic  pressure  (r2  =  0.63  —  0.95,  p  <  0.01) 
over  the  pressure  range  associated  with  breast  tumors  (0-50  mmHg)  when  measured  in 
vitro  employing  the  new  water-tank  (with  extra  acoustic  absorbers  incorporated).  This 
work  has  been  partially  published  [15]. 

Our  attempts  to  design  a  realistic  simulation  model  accounting  for  the  experimental 
results  have  been  mixed  and  further  work  is  ongoing.  These  efforts  have  been  published 
[33-34].  Software  for  analyzing  RF  data  from  the  Sonix  RP  scanner  to  produce  SHAPE 
pressure  estimates  has  been  successfully  optimized  and  several  publications  and  a  patent 
have  been  published  or  are  in  review  [32,  35-38]. 

Moreover,  in  vivo  proof-of-concept  for  SHAPE  as  a  noninvasive  monitor  of  IFP  has  been 
provided  in  5  swine  with  naturally  occurring  melanoma.  SHAPE  showed  excellent 
correlation  with  IFP  values  obtained  in  nonnal  tissues  and  in  the  tumor  using  intra- 
compartmental  pressure  measurements  with  the  wick-in-a-needle  technique  (r“  =  0.67  — 
0.96,  p  <  0.01)  [39-40].  The  work  on  in  vivo  SHAPE  measurements  of  IFP  was  selected 
for  the  final  of  the  American  Institute  of  Ultrasound  in  Medicine  (AIUM)  2011  Young 
Investigator  Award  (for  V.  G.  Halldorsdottir)  and  was  awarded  FIRST  place  out  of  90 
abstracts  submitted  for  this  competition  [39]. 

Due  to  the  delay  caused  by  the  in  vitro  experiments  and  some  personal  issues  over  the 
summer  of  2012,  the  project  was  granted  a  16  months  no  cost  extension  during  which 
task  3  was  to  be  completed.  However,  it  took  markedly  longer  than  expected  (over  6 
months)  to  arrive  at  a  stable  and  consistent  level  of  cell  growth  for  the  MDA-MB-231 
breast  cancer  cells,  which  also  produced  tumor  growth  in  vivo.  During  this  time  period 
we  developed  an  automated  AO  optimization  algorithm  to  identify  the  AO  level  eliciting 
growth  stage  subharmonic  signals  [21],  In  total,  ninety-seven  (97)  female,  nude  rats  were 
implanted  with  MDA-MB-23 1  cells  and  5 1  animals  were  successfully  studied  with  SHI 
and  SHAPE  21,  24  or  28  days  after  tumor  inoculation.  However,  the  comparisons 
between  SHAPE  and  Stryker  IFP  measurements  have  yet  to  be  completely  finalized. 
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The  SHI  images  of  tumor  vascularity  were  compared  to  immunohistochemical  markers  of 
angiogenesis  (CD31,  COX-2  and  VEGF).  The  strongest  correlation  determined  by  linear 
regression  in  this  breast  cancer  model  was  between  SHI  and  percent  area  stained  with 
VEGF  (r  =  0.53;  Table  1).  However,  we  are  in  the  process  of  re-analyzing  this  data  over 
smaller  ROI’s  to  better  account  for  tumor  heterogeneity.  Nonetheless,  this  work  on 
comparing  in  vivo  SHI  measurements  of  tumor  vascularity  in  murine  was  selected  for  the 
final  of  the  AIUM  2013  Young  Investigator  Award  (for  A.  Marshall)  as  one  out  of  11 
selected  from  a  pool  of  1 1 1  abstracts  submitted  for  this  competition  [431. 

Given  the  results  above  and  the  fact  that  our  group  has  published  twenty-seven  (27) 
abstracts  at  national  and  international  conferences  as  well  as  ten  (10)  peer-reviewed 
papers  and  one  (1)  patent  application  in  the  almost  4  years  the  original  BRCP  award  have 
been  underway,  we  submitted  an  IDEA  Expansion  proposal  to  the  BRCP.  This  proposal 
build  on  the  work  achieved  in  our  prior  BCRP  awards  by  perfonning  the  all  important 
translation  of  these  methods  towards  clinical  feasibility  by  conducting  a  first-in-humans 
clinical  trial  of  3D  SHAPE  in  women  undergoing  neoadjuvant  chemotherapy.  This 
proposal  was  funded  in  2012  and  is  currently  underway. 

In  summary,  tasks  1  and  2  have  been  completed  while  task  3  is  ongoing  (due  to  the 
parental  issues  involving  the  2  students  on  this  project). 
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Appendix  I 


The  Statement  of  Work  from  the  original  proposal: 

Objective  1 

Task  1:  Computer  modeling  and  in  vitro  experiments  (months  1-6) 

a.  Construct  an  in  vitro  experimental  pulse-echo  system  for  investigating  the  effect  of 
hydrostatic  pressure  variations  on  contrast  microbubbles  and  measuring  the 
resulting  changes  in  backscattering  (Month  1). 

b.  Design  and  modify  numerical  codes  for  a  theoretical  model  describing  the 
dynamics  of  contrast  microbubbles  under  different  pressure  conditions  (Months  1  - 
3). 

c.  Calculate  the  behavior  of  individual  contrast  microbubble  and  the  collective 
behavior  of  contrast  microbubble  populations  (Months  3  -  6). 

d.  Measure  changes  in  backscattered  fundamental,  second  and  subharmonic  signals 
for  different  contrast  agents  as  a  function  of  pressure  (Months  2  -  6). 

e.  Predict  optimal  contrast  agents  for  SHAPE  according  to  the  numerical  simulations 
(Month  6). 

f.  Select  optimal  contrast  agent(s)  for  SHI  and  SHAPE.  The  selection  will  mainly  be 
based  on  experimental  measurements  (Month  6). 

Objectives  2-3 

Task  2:  Design  and  implementation  of  SHAPE  on  a  commercial  US  scanner  (months 
7-12) 

a.  Optimize  SHI  and  SHAPE,  based  on  in  vitro  measurements  and  simulations  using 
the  actual  parameters  of  the  designated  transducers  (Months  7  -  8). 

b.  Modify  a  state-of-the-art  US  imaging  system  (the  Sonix  RP)  to  incorporate  the  SHI 
contrast  imaging  modality  and  to  perfonn  SHAPE  (Months  8  -  10). 

c.  Evaluate  the  new  imaging  modality  and  SHAPE  in  an  in  vitro  phantom  using  the 
modified  US  scanner  (Months  11  -  12). 

d.  Prepare  regulatory  review  and  obtain  approval  for  animal  studies  (Months  9  -  12). 

Objectives  3-5 

Task  3:  Animal  experiments,  data  collection  and  analysis  (months  13-  36) 

a.  Create  and  grow  breast  tumors  by  implanting  one  of  three  human  breast  cancer  cell 
lines  (SKBR3,  BT474  or  MCF-7)  into  the  mammary  fat  pad  of  athymic,  nude  rats 
(Months  13  -  34). 

b.  Calibrate  in  vivo  SHAPE  results  based  on  IFP  measurements  obtained  with  the 
intra-compartmental  pressure  monitor  in  2 1  nude  rats.  Three  groups  (one  per  cell 
line)  of  7  rats  with  breast  tumors  implanted  will  be  studied  (months  14  -  16). 
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c.  Produce  and  evaluate  the  ability  of  SHI  to  depict  nonnal  vascularity  as  well  as 
breast  tumor  angiogenesis  in  human  xenografts  implanted  in  nude  rats  compared  to 
CD31  stained  specimens  (Months  17  -  34). 

d.  Validate  the  clinical  potential  of  SHAPE  as  a  therapy  monitoring  tool  by  studying 
180  human  xenograft  breast  tumors  in  nude  rats  (42  normal  rats  and  138  after 
administration  of  a  chemotherapy  agent  paclitaxel)  and  comparing  results  to  intra- 
compartmental  pressure  measurements  (months  17  -  34). 

e.  Perfonn  statistical  analyses  and  write  final  report  (months  34  -  36). 
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